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ABSTRACT 


A method of measuring laser beam extinction in the 
atmosphere over the ocean has been devised and satis- 
factorily utilized. An attempt has been made to correlate 
the observed values with existing meteorological conditions. 
The apparatus consisted of equipment, already being utilized 
for scintillation measurements, modified by the addition 
of a source intensity monitor. Several calibration trials 
and three over-water experiments were conducted. A value 
oa 


on the order of 5x10-¢ meters was found for the 


extinction coefficient in clear weather. 
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I. PROBLEM DESCRIPTION 


A. GENERAL 

The increasing utilization of laser technology in long- 
range applications such as communications links, weapons 
designators/range finders, and high energy laser weapons 
makes it desirable to have the capability of predicting 
the effectiveness of the particular device at the termina- 
tion of the path through the atmosphere. A highly collimated 
beam of monochromatic light, such as laser light, is affected 
by many phenomena among which scintillation, beam wander, 
beam spread, and attenuation are only a few. The develop- 
ment of a procedure by which the latter, commonly called 
extinction, could be measured in a path through the atmos- 
phere directly above the ocean surface was the goal of this 
experiment. 

Atmospheric attenuation (extinction) is generally 
caused by two phenomena: absorption by the gaseous mole- 
cules, and scattering by haze, fog, rain, and other aerosols. 
The absorption problem is minimized during system design 
by selecting a lasing medium that produces light of the 


proper wavelength so as to be within one of the several 





transmission "windows" in the atmosphere. Thus, the degree 
of extinction becomes a function of the meteorological 
conditions. 

Ideally, when sufficient data has been obtained, a 
system operator would only have to observe the existing 
meteorological conditions and thus predict the effective- 
ness of his equipment. Extinction measurements over 
terrain have previously been made /Ref. 1/. However, such 
measurements over the marine boundary are virtually non- 


existent. 


B. MATHEMATICAL DESCRIPTION 
If a beam of monochromatic flux enters an attenuating 
layer, the Bouger-Lambert Law /Ref. 2/ states that the 
flux at any point within the layer is given by 
P= Rew” (1) 


where P= Flux remaining after the beam 
has traversed a path length z, 


P, = Flux entering the attenuating 
layer, 


Extinction Coefficient, 


je 
x Path length in the attenuating 
layer. 

The flux described in equation (1) above has the units 


of energy per unit time over a given area and will be 


further described. 
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Considering both absorption and scattering the extinc- 
tion coefficient 


becomes = LG (2) 


where «<< absorption coefficient, 


and x scattering coefficient. 

In this experiment, 0.6328 micrometer wavelength laser 
light was utilized and since this is within one of the 
transmission "windows" no attempt was made to distinguish 
the relative values of ~« and JY , but only to determine 
their combined value. 

The relationship (Equation 1) holds for a spherical 
wave. The detailed mathematical development, following 
references 3 and 4, presented in Appendix A, shows that a 
laser source producing pure monochromatic radiation can be 
considered to be a point source of spherical waves, pro- 
vided the observations are made in the paraxial region, at 
a distance from the source much greater than the off-axis 
position. 

Thus, considering again equation (1), and assuming that 
the effective source of the radiation is at the Gaussian 
waist of the diverged beam and radiates P, watts into a 
solid angle ™ , then the intensity at a position near the 


beam axis at a distance R is given by 


pail 





-~A-R 


T= tec» (3) 
Zoe RK 


where I = Intensity at R in Watreymeeer 
_c. R* = Area of beam at distance ee 
Consider propagation in a controlled atmosphere where 
scattering aerosols are at a minimum, so that “ = 4 , 


and a calibration factor can be determined from 


.- eas. (4) 


nm R?2 
where the subscript indicates all measurements taken during 


the controlled trial. Then from (4) 


Pe oR 
a = I,R4 hee id (5) 


2 Ae Ko MEGKo 


co t-te (6) 
R 


Taking the Natural Logarithm results in 
In (IR*) = In (1oR5 27 § - wer. (7) 
A plot of ln (TR) versus R will yield the ordinate 


intercept 
In quire ete Ry = Inc, (8) 


Le 





So in general for longer ranges where the determination 


of ue is desired 


I= -=.— (9) 
R2 
and In IR? = I1nC ~“-R (10) 
In C - In TR? 
from which: eX = ———— (11) 


Values of In S1r*8 at each value of R can thus yield 
aecalculated value of “. 

Alternatively, it is only necessary to plot In f tr2¢ 
versus R for a series of measurements at increasing ranges 
and from this determine .% as the absolute value of the 
Slope. This assumes, of course, that 4 remains constant 
over the time and distance of the measurements. 

To obtain the necessary data for the plot of In i ir2§ 
versus R the procedure briefly described below is followed: 

A horizontally fanned beam of constant intensity sweeps 
vertically across a detector located at a known distance R, 
4000 times each second. The peak value of each detected 
pulse is recorded and the pulse height distribution 
(spread out because of scintillation) is obtained. The 
average value of the pulse height is calculated and taken 


to be a measure of the average intensity I, at the distance 


R. 
13 





Pi. ~EXEERIMENTAL APPARATUS 


A. GENERAL 
Since the ultimate goal of the experiment was to measure 
the extinction coefficient over the marine boundary in con- 
junction with simultaneous meteorological measurements, the 
use of the Research Vessel ACANIA was necessary. Because 
of the limited availability of the ship's services, it was 
deemed desirable to test equipment and procedures thoroughly 
prior to actual shipboard measurements. This was accomplished 
first in the 120 meter long basement hallway of the Naval 
Postgraduate School's, Spanagel Hall. To accomplish this 
the laser source and scanning equipment (to be described 
in detail in the next section) were placed at one end of 
the hallway, and the detector and processing equipment were 
moved the length of the hallway in twenty meter increments. 
When operation of the equipment in this controlled 
environment was deemed satisfactory the apparatus was re- 
configured. The laser source and scanning equipment were 
mounted on a training and elevating mechanism in a two- 
wheeled trailer. The detecting and processing equipment 


was mounted in an instrumented bus (Fig. 1) which provided 


14 
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Naval Postgraduate ‘School Mobile Optical 
Research Laboratory. 
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Interior, Showing 18-Inch Telescope. 





FIGURE 1. INSTRUMENTED BUS and 18 <Inch Telescope 
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portability, shelter, and 110 V.A.C. power from a 3.5 K.W. 
gasoline engine driven generator. This configuration pro- 
vided the necessary support for extending the path length 

to 300 meters in 100 meter increments along O'HARE Avenue 

at the Naval Postgraduate School Annex. Operation in this 
manner was the final preliminary check prior to shipboard 
operation. For shipboard measurements the apparatus was 
again re-configured. The laser source and scanning equip- 
ment were mounted in the bus and the detecting and processing 
equipment were aboard ship. The reasoning behind this final 
re-configuration will become apparent in the following 


paragraphs. 


B. DETAILS OF EXTINCTION EQUIPMENT 

The laser source, scanning equipment, detection, and 
processing equipment utilized for extinction measurements 
were also utilized for scintillation and aperture averaging 
measurements; however, only that equipment in direct support 
of the extinction experiment will be described in detail. 

l. Bus Mounted Equipment 

Equipment described in the following paragraphs 

was mounted aboard the instrumented bus utilized by the 


optical propagation group. The bus itself was a converted 
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General Motors Corporation passenger bus. The forward 
one-third of the vehicle contained storage cabinets and 
workbench space and the drivers compartment. The after 
two-thirds of the vehicle contained electronic equipment 
racks, workbench, and an 18-inch telescope mounted on a 
platform capable of being raised clear of the floor on 
three hydraulic camper jacks extending through the floor 
to the ground. This arrangement provided a stable mount 
' for the telescope and all equipment mounted on it and was 
thus isolated from any vibration caused by personnel moving 
about in the vehicle. A portable 3.5 kilowatt generator 
provided 110 V.A.C. power to all equipment. 

a. Laser Source 

The laser utilized as a source for the extinc- 

tion measurements was a 0.6328 micrometer He-Ne laser 
manufactured by C.W. Radiation INC. The output of the laser 
was fitted with a diverging lens controlled by a micrometer 
screw adjustment. The laser, scanning equipment, and source 
intensity monitor were mounted on a common plate attached 
to the right-hand-side of the bus mounted 18-inch telescope 
(Fig. 1). This provided a means by which the laser could 
be trained, elevated, and aimed at the shipboard mounted 


detectors. 


Ie 


b. Scanning Equipment 

(1) Physical Description. Mounted directly 
in front of the laser was a mirror positioned at a 45° 
downward angle (Fig. 2). This mirror performed two func- 
tions. First, it directed the beam down onto the scanning 
mirror. Secondly, it was capable of being bent, by means 
of a micrometer screw-adjustable knife-edge placed in the 
vertical axis, so as to produce a "'FANNED" beam elongated 
in the horizontal plane. This arrangement was capable of 
producing a beam at the target ship sufficiently wide so as 
to provide uniform coverage of the detector as the ship 
moved and swung about its anchor. 

The actual scanner consisted of a small, 

1.0 cm x 1.0 cm, mirror mounted at a 45° upward angle 
[Ref. 6, p. 40/7. This mirror, attached to two piano wires 
and mounted between the poles of a strong magnet, was set 
into oscillation and driven by a power amplifier. The 
amplifier controls were adjustable so as to provide a 
maximum vertical scan angle of 14 milliradians at a fre- 
quency of approximately 2000 HZ. The result of this beam 
processing was to allow detection of maximum intensity 


independent of the ship's vertical motion. 
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(2) Mathematical Description. The following 


is a brief analysis of the scanned beam in terms of the 
vertical angular width of the beam and the scanning rate 
of the oscillating mirror. This permitted a prediction of 
the time width of the detected pulse which was then found 
to be in agreement with the measured value. To determine 
this the beam was fanned and both the unscanned and scanned 
beam sizes were measured at a convenient distance from the 
scanning mirror. Refer to Figure 3 for the details of the 
geometry; where (the distance from the scan mirror) = 
2316.5 cm, hy (the unscanned beam height at ye ) = 5.5 cn, 
and hy (the scanned beam height at f ) =P 2 Oem, 

From this ce , the unscanned beam angle 
of coverage, equals 2.37 milliradians and 26, the angle of 
coverage of the scanned beam, equals 11.22 milliradians. 
Thus, ®& = 5.61 milliradians. 


Now the pulse width “? is given by 


Pees 


/ ey (12) 


where [3 = angle subtended at the 
scanning mirror by the 
detector, 


and © = maximum angular scan rate. 
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FIGURE 3 


Geometry of Pulse Width Determination 





a ob 
Since ie = 7 for 4 >> D-~ where 
1B, equals the height of the detector. 


| Thus, 


D 
ig, ete (13) 
os 


Max 





For large values of A , the second term 


in the numerator of (13) becomes negligible thus rendering 


| 


_ the time width of the pulse nearly constant with range. 





_ : x (eS t ) 
=o ° orn GD 14 
& oe ; (14) 
or 
GO = OO Cos wT. GES) 
For maximum value of ro ; Ces wot = 
and 
: _ 16 
Frome _ IS, ) ( ) 
with =e xe DO0O At, » 
and © = 5.61 milliradians, as before; 
Thus, Sona = 70.5 radians/sec. 


For a detector of two centimeter diameter 
and yy as before, equation (13) yields a time width of the 
pulse, “ , of 45.8 microseconds. This agreed with 


oscilloscope measurements of the detected pulses. 
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Figure 4 shows the scanning equipment and 
mounting details; the entire assembly was referred to as 
the "FAN-SCAN." 

c. Source Intensity Monitor 

Periodically throughout the extinction measure- 
ments it was necessary to sample the laser output as the 
beam left the scanning mirror to ensure that no intensity 
degradation was present at the source. To facilitate this, 
the Source Intensity Monitor was devised (Fig. 4). An 
adjustable mirror was installed which could be swung into 
position to interrupt the outgoing beam and direct it back 
to the monitor detector mounted above and behind the source 
laser. The detector, a PIN-5D Photodiode was mounted in 
an adjustable mount to permit fine alignment of the sensi- 
tive surface with the sampled beam. The output of the 
detector was displayed on a Monsanto OS 226 P, oscilloscope. 

2. Shipboard Mounted Equipment 
The detection and signal processing equipment aboard 
the R.V. ACANTA included, detector, gyro-stabilized detector 
mount, PAR 113 Low-Noise Differential Preamplifier, Demodu- 
lator, Hewlett-Packard 7562A Log Voltmeter, Victoreen PIP- 


400 Multichannel Analyzer, and Hewlett-Packard HP 9810 











Laser Source and Scanning Equipment Mounted 
on 18¢Inch Telescope. 








Close-up Showing, from lef ' onitor Detector, 
Laser,Cylindrical Mirror, Oscillating Mirror, 
and Retractable Mirror. 


FIGURE 4. FAN~SCAN and Source Intensity Monitor Arrangement 
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computer and associated plotter. The following paragraphs 
describe in detail the arrangement and function of each of 
the above. 

a. Detector 

The detector utilized for the extinction measure- 
ments was a TIX L69 Silicon Avalanche Photodiode reversed 
biased to 159 volts. This provided an acceptable Signal- 
To-Noise Ratio at the extended ranges. The detector was 
mounted in a holder fitted with lens and an aperture adjust- 
able to 4.2 cm diameter. For more detailed discussion of 
the detector see /Ref. 7/. 

b. Gyro Stabilizer 

The detector was mounted on the gyro stabilized 
mount. The stabilization was necessary to maintain align- 
ment of the detector with the laser source ashore. 

The gyro was an inertially stabilized platform 
utilizing two iron flywheels mounted on each end of a 1/2 hp 
capacitor start electric motor (Fig. 5). The shore based 
laser source utilized for the extinction measurements also 
served as the beacon onto which the gyro feedback circuit 
locked. Feedback was provided by a four-quadrant position 


and intensity photodetector with associated electronics. 


> 








Gyro unit-showing top mounted Avalanche 
photodiode detector. 





Signal processing equipment; PAR113's(left); 
equipment in the rack( top to bottom) ; Demodula- 
tor, Log Voltmeter/Amplifier, Oscilloscope,and 
Multi-channel Analyzer. 


Figure 5. Signal Detection and Processing Equipment 
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A small aiming telescope and "JOY STICK" control 
were mounted on the left side of the — to aid the gyro 
operator in aiming the control detector at the laser source. 
The gyro platform was positioned on the upper deck of ACANIA 
and was sheltered from the elements by an aluminum tool 
shed open on one side. The signal processing equipment 

| described below was bench mounted in the "DRY LAB" on the 
ean deck. (Fig. 5). 
c. Preamplifier 
The output of the detector was fed through an 
| E.T. follower (which served the function of impedance 
matching) to the PAR 113 Low-Noise Preamplifier. This 
provided high gain, 10 to 104, low noise amplification, of 
the received signal. Adjustable High and Low Frequency 
Rolloffs allowed selection of the bandwidth. 
d. Demodulator 

The modulated signal was next passed to the 
Demodulator. This device sensed the scanning beam pulse, 
followed the pulse to its maximum and held that instanta- 
neous signal level in a flat-topped pulse until the next 
maximum scanned across the detector at which time it held 
that maximum level, and so on. The output of this circuit 


was thus a series of step functions. 


Za 


e. Log Voltmeter 
The Hewlett-Packard 7562-A Log Voltmeter was 
the next instrument in the processing step. At this point 
the Logarithm of the signal was obtained. 
f. Multichannel Analyzer 
The Victoreen, PIP-400, Multichannel Analyzer 
collected, stored for further readout, and displayed the 
sampled Log Voltage Pulses. Each of the available 200 
channels corresponded to a different pulse height and the 
number of pulses of any particular height received during 


the sampled period were stored in the respective channels. 





g. Computer and Plotter 

The output of the Multichannel Analyzer was 
fed into a Hewlett-Packard 9810A calculator. The calculator 
provided controlling functions to the associated plotter. 
The calculator sampled each channel of the Multichannel 
Analyzer, stored the information and plotted the pulse 
height distribution. The calculator then computed a least- 
squares polynomial fit of a Gaussian distributed curve to 
the data points. Using the position of this curve with 
respect to read-in calibration points and the Standard 
Deviation of the curve, the average intensity of the re- 


ceived signal was computed. See Appendix B for a detailed 
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analysis of the procedure for obtaining the average pulse 
height from the log voltage distribution. For a detailed 
description of the computer program see Ref. 8. Figure 6 
is a sample printout from the system. 
h. Humidity Measuring Equipment 
The Hygrodynamics Digital I hygrometer indicator 

and the Dunmore-type Lithium Chloride sensor were used to 
measure the Relative Humidity. The equipment operates on 


the basis of resistance change within an electrolytic 





solution generating a reference voltage variance which is 
) proportional to the relative humidity change. For further 


| @erails see Ref. 9. 
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tt. “SHE EXPERIMENT 


A. EXPERIMENTAL GOALS 
The original goals of the experiment were threefold. 














|First, the main goal was the development of the equipment 





and procedures to measure extinction and to demonstrate 
that extinction could be measured in the ocean environment. 
Secondly, it was desirable to correlate the observed ex- 
tinction with the measured quantities of airbomrme particu- 
late matter and water vapor measured as humidity. Thirdly, 
it was desired to measure the extinction caused by fog and 
Pan. This last goal was not achieved since, during the 
scheduled shipboard trials, the weather was persistently 


clear. 


B. EXPERIMENTAL PROCEDURE 
1. 22-23 January 1976 
On these dates the first equipment trials were 
performed in the basement hallway of Spanagel Hall, Naval 
Postgraduate School. The laser source and monitoring 
equipment were mounted on a portable equipment rack and 
positioned at distances of 15.8 m, 31.0 m, and 46.3 m from 


the receiving equipment which was kept stationary. The 


Sil 


purpose of these trials was to develop technique and 


initial equipment settings. At that time the signal 


' detector was a PIN Photodiode which, for reasons to be 


explained in Section III.B.5, was replaced by the Avalanche 


Photodiode previously described. In addition to the devel- 


opment of techniques and equipment settings, these trials 


established the Laser Source Strength Monitor reference 


reading at 31 mv. 


2. 26 January 1976 


On this date the previous trials were repeated with 


similar results. In addition, the effects of ''Fan-Scan" 








amplifier settings and PAR amplifier settings were measured. 


It was determined that "Fan-Scan" amplitude had a definite 


effect on the monitor pulse height and an amplifier setting 


of 1.2 was chosen as the best setting since the response 


at settings surrounding this point was relatively constant. 


The PAR settings arrived at were: (1) Low Frequency Rolloff 


100 HZ, (2) High Frequency Rolloff 300 KHZ. These settings 


provided the cleanest signal response with minimum pulse 


shape distortion. 


3. 2 February 1976 


On this date, the first concrete attempt was made 


to determine the equipment response to transmission through 


a2 





an atmosphere containing essentially no water droplets. 


‘The trial was again made in the basement hallway at trans- 


: 


and 122.5 m. All equipment settings were maintained as 


Meeesion distances of 15.8 m, 31.0 m, 46.2 m, 61.5 m, 92.0 m 


previously described; however, in this trial the laser 





source was left stationary and the detector and signal 
processing equipment were positioned at the prescribed 
ranges. This procedure was utilized since it facilitated 
more rapid detector alignment. The results of this trial 
i plotted in Figure 7. The observed value of in J 1R2 
determined here was to be utilized as In {C§ in the cal- 
culation of the extinction coefficient as a result of the 
data to be taken during the shipboard trials of the week 
o£ 29 March 1976. 
4. 18 March 1976 

It was decided that, prior to the 29 March Acania 
‘trip, it would be advantageous to have some measurements 
at longer distances than could be achieved in the basement 
hallway. For the calibration trial on this date the laser 
source and monitoring equipment was mounted in the small 
two-wheeled trailer, and positioned at the western end of 
O'Hare Avenue (NAVY ANNEX) in close proximity to the fuel 


station so that 115 volt A.C. power could be readily 


33 


i ‘s 
»e% St 
vy reer 


- © 8a Be 








O0€ 





9161 Azenzqegq 7 JOZJ Y “SA (-UL)UT JO udeAD =f TANIA 


SHJ1IW NIJONYY 
002 00L 0 


9/61 Aasensqsy Z 
juawmaseg jjeH jsheueds 9 
Byeg wuomesgiyje) 





34 





stat neliondtis2 


jvenssed ite lageang? 
vie? vrveruisl S 


| 


| 


obtained from the distribution panel located there. (See 





Figure 8). The detector and signal processing equipment 
was located in the instrumentation bus, which was moved 
and positioned at ranges of 91.0 m, 206.0 m, and 336.0 m 
‘from the source. (See also Figure 8). The results of this 
trial are plotted in Figure 9. These results together with 
ithe results of the 2 February trials were to be utilized 
with the data collected aboard Acania during the week of 

29 March. 


| 


| 5. Week Beginning 29 March 1976 
| During this week, several extinction measurements 
were recorded aboard R. V. Acania. The majority were taken 
with the ship at anchor at a range of approximately 0.8 
nautical miles off Point Pinos. Each day the ship would 
arrive on station and anchor using radar ranges. When the 
ship settled and rode to the anchor, the exact range from 
bus to ship was determined by triangulation using surveyor's 
transits. It was during this trial, on 30 March 1976, when 
an attempt was made to obtain extinction measurements at 
consistently increasing ranges, that it was determined that 
to do so would require the use of the Avalanche Detector; 
Since, at ranges greater than two kilometers, the Signal- 


to-Noise Ratio obtained with the PIN Photodiode was 
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insufficient. From this determination, it was decided that 
a "calibration" trial, similar to the 18 March trial, would 
have to be made using the Avalanche Photodiode. Since this 
detector produces an extremely large signal resulting in 
‘Saturation at close ranges, it was decided to again work 
hat the Navy Annex site rather than the Spanagel Hall basement. 
6-19-20 April 1976 

On these dates the "calibration" trials were con- 
ducted utilizing the Avalanche Photodiode. The equipment 

bak again configured as described previously in Paragraph 4. 
The source strength monitor at this time also produced a 
reading of 31 mv indicating that the laser output was re- 
maining constant. The results of the trials described here 
are plotted in Figure 10. Again, the observed value of 

In 31R*{ determined here was to be utilized as ln $C 

in conjunction with the data to be collected during the 

final shipboard trials to be conducted the following week. 
Also, during this calibration trial, anticipating the need 

to use a large aperture opening at the longer ranges, a 

trial was conducted to determine the effect on the response 

of the detector as a function of aperture size. This was 


conducted at the longest available range (336.0 m) since 


opening the aperture greater than 0.4 cm diameter at the 
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closer ranges resulted in saturation. The results of this 
‘trial are plotted in Figure 11. As indicated the response 
is linear to an aperture radius of 1.1 cm, at which point 
saturation set in. 
7. 28-29 April 1976 

On these dates, the final data collection reported 
on in this thesis was completed. Simultaneous meteorological 
data was collected at the detector site aboard ship. The 
laser source was mounted on the 18 inch telescope in the 
| 
bus. The bus was positioned at the water's edge on Point 
hee. The ship was positioned north of this location at 
various ranges from 839 meters to 6/748 meters. Three 


factors arose that affected the effort to duplicate the 


exact equipment conditions of the calibration trials. First, 
| 


because of the heavy electrical load on the 3.5 kw generator 





resulting in a low line voltage, the ''Fan-Scan" amplifier 
setting had to be increased to 2.5 to maintain oscillation 
of the scanning mirror. Secondly, the source intensity 
monitor indicated some laser degradation as a result of 
the low voltage condition. The maximum monitor reading 
was 22.0 mv when the trial first began, then within 30 
minutes of operation the monitor dropped to 18.5 mv and 


remained there throughout the remainder of the trials. 
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This was down from a constant 31.0 mv obtained during the 


calibration trials. Thirdly, because of the long ranges 


involved, it was decided that the best procedure would be 





to set the aperture size at 4.0 cm diameter and keep that 





constant setting throughout all remaining trials. 





The equipment settings were as follows: 
| PAR Amplifier: Low Frequency Roll Off, 100 Hz., 
High Frequency Roll Off, 300 kHz; 
"Fan-Scan" : Amplifies:, 2.5, 
Cylindrical mirror micrometer, 7; 
Laser divergence micrometer, 5; 
Detector Aperture, 4.0 cm diameter. 
The general experimental procedure for these two 
days was to maneuver the ship into position north of the 
bus at the desired range. When on station the gyro stabi- 
lized detector would be aimed at the laser source. When a 
usable signal was obtained the data collection would 
commence. The Multichannel Analyzer would be allowed to 
collect a sufficient number of counts so as to produce a 
nearly Gaussian distribution. This usually took from one 
to three minutes depending upon distance from the source. 
Since 1t was nearly impossible to hold the ship at the 


same range throughout data collection, approximately half 
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jway through the collection process the ship's range from 
idee source would be determined by triangulation with the 
eror' 2 transits. As soon as the data collection was 
completed the ship would transit to the next desired 
position and the process would be repeated. While the ship 
was moving to the new position the data was fed into the 

HP 9810 Computer, processed, and plotted on the associated 
plotter; thus freeing the Multichannel Analyzer to collect 
data at the next position. At the closer ranges the gyro 
would "lock-on" the laser signal; however, at ranges greater 
than about 3000 meters the signal strength was not sufficient 
for "lock-on" thus the gyro operator controlled the aiming 
lof the detector manually with the '"Joy-Stick'’ control. The 
results of the data collected and the correction factors 


resulting from the equipment changes previously described 


are reported in Section IV, RESULTS AND CONCLUSIONS. 
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IV. RESULTS AND CONCLUSIONS 


A. GENERAL 

Since the major goal of the experiment, as stated 
earlier, was the determination of the extinction coefficient 
within the marine boundary layer, the actual data gathering 
lewoived about the schedule of the Research Vessel Acania. 
—s of this and the discovery on 30 March of the long 
range limitations of the PIN Diode detector, the data 
breaks down logically into two distinct experimental evolu- 
‘tions: (1) the 30 March experiment and trials leading up 
‘to it, and (2) the 28-29 April experiment and the calibra- 
‘tion trials leading to it. 
| 

To facilitate the determination of the extinction co- 
efficient, xc , considering the paucity of long-range data, 
the Natural Logarithm of the quantity (IR) was plotted 
versus the range, R, for each data point obtained. A 
straight line was drawn through these data points utilizing 
the ordinate intercept and slope as determined from the 
Principle of Least Squares curve fitting /Ref. 10/. The 
absolute value of the slope of this curve is equal to 4c ; 


the extinction coefficient. 


44 





The calibration trials served to provide short-range 
data and thus determine, to a large degree, the ordinate 
intercept of the plot; whereas, the long-range data obtained 
in the desired ocean environment, from Point Pinos to Acania, 
provided the mikaeemel on for the slope and consequently ,<« . 


| 
Prior to plotting the data, the following corrections 


| 
| 
were made where necessary: 


1. Source Intensity Correction 
Since the majority of the data was taken with a 
‘laser source intensity monitor reading of 31.0 mv, the 
eat height data taken while the laser intensity monitor 
reading differed from this was corrected either up or down 
according to the ratio of the two readings. 
2. Aperture Size Correction 
As mentioned earlier, it was necessary to utilize 
a small aperture at short ranges and a larger aperture at 
long ranges. Therefore, the data collected at the longer 
ranges was reduced by the ratio of the squares of the 


aperture radii. 


B. THE 30 MARCH EXPERIMENT 
The 30 March 1976 experiment and the calibration trials 
preceeding it served to establish the feasibility of making 


measurements at increasing ranges as well as providing the 


45 





initial measurement of the extinction coefficient in the 
marine environment. The data is tabulated in Appendix C. 
Figure 12 is the plotted data including the least squares 


| 
fitted straight line. The weather was clear and the simul- 


taneous meteorological observations provided an average 


relative humidity reading of 0.911 over the time of the 


shipboard measurements. At was 6.2 x 1074 meters — 


c, THE 28-29 APRIL EXPERIMENT 

The experiments conducted on 28 and 29 April 1976, 
together with the preceeding calibration trials, succeeded 
in providing data at even longer ranges. The data taken 
on these two dates was considered together for calculation 
purposes because of the similar weather conditions and 
equipment operation. The data is tabulated in Appendix C, 
Figure 13 is the plotted data including the least squares 
fitted straight line. The weather was again clear and the 
relative humidity on 28 April was 0.776, while on 29 April 
it was 0.774 giving an average relative humidity reading 
of 0.775 over the time of the measurements. << was 5. ome 


1074 meters ~-, 


46 


000% 





9/61 YOIRW OF AOF WY “SA (,UI)UL JO ydexzy 


SHYJLIW NI JONWY 
OO0€ 0002 


:eyeg uolesquyeg 


@® :283e89 yaseW OF 


» OLX 29— = adojs 


9/61 YI4EW OF 
VINWIW 03 SOUld Yd 


él AMNOLA 


0001 


47 





0009 


9/61 Tt4idy 67-872 FOF U “SA (7UI)UT FO ydeszp 


SHJLIW NI JONVY 
000s O00F OOOE 


¢:8389 woyesqgujeg 
Ww :e3eg jlady 672 
© : eyeg judy gz 


Ol XY gE- = adojs 


Og 


‘2 


9/61 1!44¥ 62-82 
VINYVIOY 93 Souldg ‘jd 


Cl dAaNoLTs 


0002 


0001 


48 


. a ; , a. —_ 
-« " i: : 
ris me .- = —_ ee — a -_ 
re eh 7 " : Oo 
7 _ es 
apn. “2 SF) **. 
- “— : i 
- a ie _ 
- 7 7. 3 7 : 
: : - a 2 
-_ 
- 
7 
_ > 
¢ a 
7 
5 7 rs oe 
See 7 : . 
FP se: _ a 
Tye 4 - 7 ] 7 7 -_ ; 7 
- ‘bgt _ 
> 7 
- a . a oe ° ¥ a 
yaw. Sh 7 
-O ite ge G} : a 
OO _ - 
_ ma 7 ore oo 7 ah _ - _ 
— a a a a - a OO 7 : - 7 
a re ee one — - ee - a 
7 _) ~~ =-ix a Lf 2. .* “ = os 





a 
—— 











— 





— — 
—_ + _ 
- ae 
7 S 
a - a 
_ 
> 
7 
- 
7 
> - 
-_ 
7 
a 
_ a _ 
_ - _ : _ 
: — _ - ; - 
; _ 
= ¢ v ~ 
aS a —— ian a 





D. CONCLUSIONS 

It seems possible to draw three main conclusions from 
these results. 

First, it is experimentally feasible to measure laser 
beam extinction at Sea over the ocean-air interface. 

Secondly, the extinction coefficient for 0.6328 micro- 
meter laser light in clear weather, through the marine 
environment is of the order of 5 x 10~* meters ~!. 

Thirdly, humidity may be considered as a possible para- 
meter with which to predict laser beam extinction in clear 
weather. 

The latter two items, of course, require further 


investigation. 
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V. DISCUSSION OF ERRORS 


The considerations and assumptions discussed below were 
utilized in an attempt to understand how the uncertainties 
in the measured quantities affect the uncertainty in the 


determination of “ec - 


one —L2n 
Tie Cu 


i 
as (17) 


the more useful form 


hy ECR: LC var, (18) 


was used. Ehen 4 was determined to be the absolute value 


of the slope of the plot of ln IR* versus R. From (18) 


ae = 


n — nr -aknR 


From the theory of propagation of errors /Ref. 10/; 


(An) ” = (2) (ae)” + fae) ‘(se (34) @ R). (20) 


‘From equation (19), 


dn . | 7 en a 
eee 2 sr 8 | Re 
and DH mae. ( AtuR . Thus, the expression for 
IR Rr 
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the relative error in 7 becomes 


2. 


ac)" (AE) (at a) (28) 7 * 
5 a ae < ) (a) 


au 
Ze _<t # 


where the uncertainties in C, I, and R are assumed to be 
independent. For the = determination, two types of range 
error were assumed to exist; one because of ship drift during 


the time data was being collected, and the second because of 


surveying uncertainties. Therefore, 4% becomes 





gees (22) 


) 


where (29) is the relative error introduced by assuming a 


0.05R 
6000 


assuming a 300 meter error at 6000 meters and a lesser 





20 meter ship drift, and the ( ) term arises from 


error as range decreases. Furthermore, it is assumed that 


cs = ss = 0.05, and a value of 5.0 x ge is utilized 
for 4c : 


Then, utilizing equation (21), the above assumed values, 
and ranges from 1000 to 6000 meters in 1000 meter increments, 
a plot of the relative error versus range was constructed 
ce. 14). 

In addition to the above investigation, the least squares 


calculations used in the curve fitting process were used to 


eye 
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determine the variance in the slope of the plotted data. 
The Standard Deviations of the data for 28 April and 29 
April were found through separate calculations; the 
variance of each day's results were then computed. 

This procedure indicated a relative error in of 


16%. 


a) 





VI. SUGGESTIONS FOR FURTHER INVESTIGATION 


The following paragraphs are intended to present several 
ideas and modifications that, if implemented, would allow a 
more accurate determination of the extinction coefficient 
under the conditions experienced during this type of experi- 
ment. Furthermore, it would be possible to investigate 
the propagation of laser produced radiation at other wave- 
lengths. 

A new Scanning device, recently procured, enone serve 
to stabilize the '"Fan-Scan" output under varying line volt- 
age conditions. 

A second demodulator unit, presently under construction, 
will provide an input to a digital voltmeter, thus facili- 
tating a more accurate and rapid determination of the laser 
source intensity. Additionally, a method of monitoring 
the signal detection equipment performance would prove 
valuable. 

The use of a Laser Rangefinder should be considered. 
This would provide more accurate and rapid source-to- 
detector range determination. 

Ship to ship, or ship to some other waterborne platforn, 


operations should be considered. Since the shoreline, across 
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which the propagation path passes in the present experi- 
mental setup, is rock strewn and relatively unprotected, 
waves breaking on these rocks produce a significant amount 
of airborne water particles and vapor. This being directly 
in front of the source can only act to produce a fallacious 
result in the determination of extinction. An experiment 
conducted entirely over water would eliminate this problem. 
Further attempts should be made to measure 0.6328 micro- 
meter radiation extinction in fog, mist, and rain. Also, 
an effort should be made to meaSure the extinction coeffi- 
cient of other wavelength laser radiation; specifically, 


CO4 produced, 10.6 micrometer radiation. 
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APPENDIX A 


SPHERICAL WAVE SOLUTION TO THE WAVE EQUATION 


The following development will show that a laser, which 
produces pure monochromatic radiation, can be considered to 
be a point source of spherical waves as long as the region 
of observation of these waves is restricted to a region 
near the beam axis and at a range from the source much 
greater than the distance off the axis. 

Consider now a laser source in which the intensity 
distribution at planes normal to the propagation is 
Gaussian. Starting with the Maxwell Equations in a homo- 
geneous charge-free medium, so that Gie-o , the 


characteristics of the Gaussian Beam are developed /Ref. 3/. 


Then Ven = ¢ == +a ; (A=1)) 

Where ae = oF a 

Also Vx —-—/h (A-2) 
~~ ae 


Taking the curl of (A-2) and substituting (A-1) for 


Ta results in 
Vx (Wak)? sx slesH/, - ah Es E26 roe eh (A-3) 
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— ae 
and vax (vx E) = THE 3 2 rug c& ) (A-4) 
dt dt 
Qiae = 
ae ¢ E <- =O , (A-5) 








where VxVYxe = V(v-6)-VEéE, 
and V-o@ =C. 


Assuming that the field quantities vary as £(%, oe 2,0) = 
> ‘cot 
ie >E% 4,2) ae” iC that is pure monochromatic radiation 


such as from a laser - the Wave Equation (A-5) becomes 
-» — . =~ —> 
VE +R (M)DF 29 , (A-6) 


where ps wwe §1-4E 


and (A-6) allows for the possible dependence of k on 
position r. At this point the development is limited to 
the case in which k2(P) is given by 

Ke ( F) =k? - kor” (A-7) 


where Ko is some constant and k = Ze , where gn is the 


wavelength of a TEM (Transverse Electromagnetic) plane 
wave propagating in the medium. Considering the propaga- 


tion of a nearly plane wave in which the flow of energy is 


5/ 





predominantly along the # axis and assuming a solution to 
(A-6) whose transverse dependence is on r = |x 4 y? only. 


then vy? in (A-6) can be replaced by 











a Ae ar a eee) 
V VY. + ae oye - ror > (A-8) 
taking E as 
_ —ikz2 
ZL = WV (x, A) z) a * ) (A-9) 


equation (A-6) becomes 


5 z - S - ik 2 
y*blay,2)e** b Play,2)o~*_ kkar* Vix y,2)e" =0, (A-10) 
and operating on (A-10) with (A-8) 


Tt -2iky'-bh ry =, (ALD 


) 
where Y = IY , and the assumption is made 


G2 


that the transverse variation is slow enough that 
) u Z 
RY >> P << kl. 


Taking Y in the form of 
: l Z 
p= expf-i [Pa +tg@cls . a 
substituting into (A-11), and using (A-8) gives 
Oe 30-26 -0'-9f Peer ea 2 (A-13) 
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If (A-13) is to hold for all r, the coefficients of the 


different powers of r must be equal to zero. Thus 
2 b) 
QO” ko kia 0, (A-14) 


and p’ = = = CAS ) 


Now according to (A-7), if the medium is homogeneous, k, 


can be set equal to zero and (A-14) becomes 


Q + kQ’ = 0. (A-16) 
Now if the function S(z) is introduced by the relation 
o) 
6, B= hx Ss”? (A-17) 


It follows from (A-16) that S" = 0 and Ss =a also 
S= az + b so that from (A-17) 

Q(Z) = k ase - (A-18) 
Following Kogelnik, /Ref. 4/, it will be advantageous to 


define a parameter q (z) where 


= R_ = Bhs 
gz) - O12) 2 Oz) . (A-19) 
Now (A-18) may be rewritten as 
iia (A-20) 


Equation (A-15) may now be rewritten as 


’ 


P = 


-# Poa (A-21) 
& 2+ Zo 2 








and 





Pig = iL(irts ‘ (A-22) 
Po 
where the arbitrary constant of integration is chosen as 


zero, which amounts to fixing the phase of the final solu- 


tion at zero. Substituting (A-20) and (A-22) into (A-12) 


pnene§ni[-ibe (uf) + era 1S, wm 


Taking the arbitrary constant of integration q, to be purely 
imaginary which leads to waves whose energy density is con- 
fined to a region near the Z axis, and re-expressing it as 


a new constant We, 





2,2 ATW, (A-24) 
an 
where 
a: le (A-25) 


Using (A-24) and the fact that 
ay, b 
L, (arib) = af Arepe #1 Can (4) 


the first term of (A-23) becomes 
/ AZ 
en ee) aa ae (A-26) 
exe|-J, (-aRb))* (eRe ae ta Coan 
ai 
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Substituting (A-24) into the second term of (A-23) and 


separating into real and imaginary parts 


eae ae go _ 
exp] aera) | Cla [ir (BE J] 2112 J] a 


Defining the following parameters 





td (2) = whe i é ge.) / j (A-28) 


R= a2 ue) we (A-29) 





¢ 





ad (25 2 ; (A-30) 


TY to 


Combining (A-26) and (A-27) into (A-23), and utilizing the 


— — ike 
fact that E(% 2) = 0 (% 4, z)e results in 


as ae: k 
Etg,2) ~ 2, exp [-4 (ha-®- "(Ste * Fe) [OD 


Consider now the form of a spherical wave emitted by a 


point radiator placed at Z = 0 /Ref. 3/. It is given by 
EE «x +e (-ck Vx" ty2°) 
[ge J 


= tke-ik x+y” : 
X~ 2 2xp ( cS 44 ), (A-32) 
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2 2 


provided that X* + Y2<<Z*, the paraxial limitation, is 
observed. 

The last form of equation (A-32) follows where Z is 
equal to R, the radius of curvature of the spherical wave. 


Since r? = x2 + a equation (A-32) becomes 


= = 7 £xpP (-ikz ~ re e! 2 (A-33) 


Now, comparing the form of equation (A-31) to that of 
(A-33), it can be seen that (A-31) is in the form of a 
spherical wave, provided the observer is in the paraxial 
region at a distance much greater than the off-axis 


distance, and also provided that 22 >? i : 
TP ad, 
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APPENDIX B 
PROCEDURE FOR DETERMINING AVERAGE VOLTAGE, t7 : 
FROM LOG VOLTAGE DISTRIBUTION 
The following development will serve to illustrate the 
validity of taking the peak value of the nearly Gaussian 
Log voltage distribution, as modified by the standard devia- 
tion of the best fit Gaussian curve, to be a measure of the 


= 


average observed voltage V . 


Let 
ys Kha (B-1) 
where 
Bk = channel number in multichannel analyzer, 
V = voltage into the log voltage circuit, 
4 = Ded, (UF and Vi known calibration 
fp, ve voltages) , 
V; 
and A : Ui. 
‘ 2 #5 eai/K 
Now, define 
x= Oy (B-2) 
Po chat Ae 
x= de Us (B-3) 
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and Ah ~ ee 


we (B-4) 
From (B-3) and (B-4), 
uU= VU; a” (B-5) 
Now 
Purydyr = pllx (B-6) 
where 


{ Peele = { pliaelx 


(B-7) 
Thus, the average voltage 1 can be expressed as 
>. 
= 6 P(r) Av ) (B-8) 
or using (B-7), 
- = \ Pwredyv, 
= ) eae 
from (B-6), 
yo 
. { pwela, (B-10) 
Substituting (B-5) into (B-10) for U results in 
- 2. & 
aL Sv; ec. pliaylhx, (Beas) 
Now assume 
~ (X-Ke) 
PY=He 207 , | (B-12) 
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and replacing p(x) in (B-11) results in 


ome | ae Mee fz, (B-13) 
or 
a ._ 1 G-%@ =) oe % 
ne = 4 et FO He 
v C ¥ Lk, (B-14) 
Expanding the exponential: 
) 7 
2 ar a 
2 = 20% MTG Ors) aoe 
X-Ke) -29% XK = r & =A 2) (8-15) 





LORe Tx * 


Replacing the exponent in (B-14) by (B-15) results in 
[x-(%+G I) m+ BS 
= —e po” 
1- = ey. a Um ~—Z AX, 


x 4 
Np HE 


Simce 2, * is independent of X , it can be removed 


= (tpt zy)" 
2072 olf 


(B-16) 


from under the integral, and recognizing that Ate 


is equal to C, from (B-7) and (B-12), equation (B-16) 
becomes 
out 
cs Xp te aa 
i = Qare ; (B-17) 


In actual operation, p(x) was the distribution collected by 


the multichannel analyzer. In calculating the best Gaussian 





curve fitted to the data distribution, the computer 
generated Xp, O%, and subsequently 1 . The validity 

of this result depends solely on how well the computed 
Gaussian actually fits the experimental distribution. This 
was indicated graphically by the curve drawn onto the data 
(see Figure 6). 


Reference 8 describes the computer program in detail. 
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APPENDIX C 


DATA 
SPANAGEL HALL BASEMENT 2 February 1976 
RANGE (meters) In(IR2) 
31 6.89 
62 71s 
100 ys. 
WAZ 7 ou 
ISK, TL VAS 
200 Tals 
250 6.95 
DLIES: 6.80 
NAVY ANNEX 18 March 1976 
91 TR23 
91 7.26 
181 7.01 
181 Sysioks. 
181 6.99 
S35 7.04 
B85 2206 


67 





NAVY ANNEX 19-20 April 1976 


RANGE (meters) In(IR2) 
91. Sab9 
on Saz2 
91 8.25 
91 8.15 

206 8.14 
206 8.30 
3316 7.99 
336 8.34 
336 8.14 
PT. PINOS to ACANTIA 30 March 1976 

RANGE (Meters) In(IR2) 

1604 Dez 

2384 By laly, 

S925 baby, 

3915 4.74 
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PT. PINOS to ACANIA 28 April 1976 


RANGE (meters) 1n(IR*) 
279 Ve 2 
2195 (55 
3288 7.42 
4420 6.60 
4535 once 
5302 5S 
5425 509 
6748 Seo 

PT. PINOS to ACANIA 29 April 1976 

RANGE (meters) In(IR2) 
839 Ts 
2120 Fe fils 
SV 5 I ase 
*4156 eS ye: 
5300 G22 
6283 See 


*The following is an example of the calculations and 


corrections applied to the value of U , obtained from 
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the pulse height distribution, to arrive at the value of 
In (IR2) tabulated in this Appendix. The particular data 
point used in this example was chosen, since all correc- 


tions possible were applied in the determination of In (IR) . 


6 Average voltage from pulse height distribution 
x 


Pre-amplifier Gain Setting 


in 3 1r2$ 


Ik 











31 mv 0.04 2 
Monitor ~ (Aperture Radius)? . (Range) *f 
Reading 
0.538 31 0.04 
=a ; Co oes aa ee (4156) 74 


in f 1557.1244 6 


= 7.35 
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